Towards a Relativistic Description of Exotic Meson Decays 
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This work analyses hadronic decays of exotic mesons, with a focus on the hghtest one, the 
J^'-' = 1^+ TTi, in a fully relativistic formalism, and makes comparisons with non-relativistic results. 
We also discuss Coulomb gauge decays of normal mesons that proceed through their hybrid compo- 
nents. The relativistic spin wave functions of mesons and hybrids are constructed based on unitary 
representations of the Lorentz group. The radial wave functions are obtained from phenomenolog- 
ical considerations of the mass operator. Fully relativistic results (with Wigner rotations) differ 
significantly from non-relativistic ones. We also find that the decay channels tti — > yrfei, tt/i, KKi 
are favored, in agreement with results obtained using other models. 
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I. INTRODUCTION 

Strong hadronic decays have been a subject ofphe- 
nomenological investigations for several years Q, |^ |M 1^ 
0, 1^ . The quark model description of strong decays is 
based on the assumption that decays originate from cre- 
ation of a quark-antiquark (QQ) pair in the gluonic field 
of the decaying meson. The produced QQ pair subse- 
quently recombines with the spectator constituents and 
hadronizes into the final state decay products. Such a 
picture is consistent with the experimental observation 
that most hadronic decays involve a minimal number of 
final state particles in the final states; this requires a 
small number of internal transitions at the quark level. 
Furthermore, the absence of complicated multi-particle 
hadronic decays is also consistent with a minimal analyt- 
icity assumption for the scattering amplitude. The ab- 
sence of complicated multi-particle cuts in the s-channel 
provides justification for a simple resonance Regge pole 
parametrization of the t-channel amplitude over a wide 
s-channel energy range. 

The mechanism of QQ pair production is by itself 
a complicated phenomenon, which can in principle be 
studied through lattice gauge simulations. In the strong 
coupling and/or non-relativistic limits, pair production 
was shown to be similar to electron-positron production 
in a strong uniform electric field (the Schwinger mech- 
anism) 0. In this case the produced QQ pair carries 
the quantum numbers of the vacuum, i.e., unit spin cou- 
pled to unit relative orbital angular momentum. A phe- 
nomenological hadron decay model based on such a QQ 
production mechanism is referred to as the ^Pq decay 
model (the spectroscopic notation refers to the quantum 
numbers of the produced QQ pair); this model has been 
extensively used in phcnomenological studies of meson 
and baryon decays. 

Within a non-relativistic or constituent quark model 
Born-Oppeinhemer approach, it is assumed that forma- 
tion of gluonic field distributions decouples from the dy- 
namics of the slowly moving constituent quarks. Con- 



sequently the decay and formation of final state hadrons 
can be studied within a non-relativistic quantum mechan- 
ical framework jSj. For light quarks the non-relativistic 
approximation can be justified from the observation that 
dynamical chiral symmetry breaking leads to massive 
constituent quarks and transverse gluon excitations with 
a mass gap on the order of 1 GeV [|i . This heavy, 
effective mass of gluonic excitations results in weak mix- 
ing between the valence and multi-particle Fock sectors 
and suggests the validity of an approach to the decay 
process in which the pair production interaction is used 
only once. 

This is also consistent with the characteristics of the 
experimental data on decays discussed earlier. Since fully 
dynamical lattice simulations of the light hadron reso- 
nances are not yet available, such a phcnomenological 
approach seems to be a reasonable starting point towards 
a description of the strong decays of light hadrons. 

In this paper we discuss some the remaining issues per- 
taining to models based on the ideas presented above. 
The first issue is the question of relativistic effects. Even 
though a simple non-relativistic description appears to 
be quite successful in predicting decay widths of reso- 
nance with masses as large as 2-3 GeV, the presence of 
light quarks moving with average velocities a substan- 
tial fraction of the speed of light naturally raises the 
question of the validity of the non-relativistic reduction. 
In the case of form factors, it has been demonstrated 
that relativistic effects in the quark model are in gen- 
eral quite large 0, 0, 0, 0| . A study of relativistic 
effects in light-cone quantization has also recently been 
performed .15j. In the absence of a fully consistent dy- 
namical model, one could argue that these effects might 
somehow be mocked up by effective parameters, never- 
theless if one seeks a more " universal" constituent quark 
model of hadronic properties, it is essential to address 
the role of relativistic corrections in decays. 

Secondly, as discussed above since the ^Pq decay model 
can be related to the characteristics of a Wilson loop, in 
a consistent description of a decay process one should 
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consider hadrons including flux tube degrees of freedom 
and the effects of flux tube breaking [Ta IT^ . The ma- 
jority of models of normal hadron decays do not include 
such effects. In those models the strong decay ampli- 
tude is determined by the quark model constituent wave 
function, multiplied by a form factor representing pair 
production that is independent of the quark distribution 
in the parent hadron. This would not be the case for a 
general string breaking mechanism. 

A simple phenomenological picture of hadrons and 
their decays in terms of quantum mechanical wave func- 
tions emerges naturally in a flxed gauge approach. For 
example, in the Coulomb gauge the precursor of flux tube 
dynamics, including string breaking, originates from the 
non-abelian Coulomb potential, which also determines 
the quark wave functions The string couples 

to the QQ pair via transverse gluon emission and absorp- 
tion in the standard way. This coupling carries '^Si quan- 
tum numbers. It is interesting to investigate whether this 
coupling, combined with the flux tube dynamics, is con- 
sistent with the phenomenological ^Pq picture discussed 
above. 

If one were to attempt a description of decays based 
on Coulomb gauge quantization, it is necessary to ad- 
dress the role of the hybrid quark-antiquark-gluon con- 
figurations, since these appear as intermediate states in 
the decay process as illustrated in Fig. ^ If such hybrid 
states also exist as asymptotic states they would provide 
an invaluable insight to the dynamics of confined glu- 
ons. In recent years, evidence has been presented that 
such states do indeed exist, in particular in exotic chan- 
nels that do not mix with the QQ sector 21, 22, 23J. 
Recently the non-relativistic ^Pq decay model has been 
extended to study hybrid meson decay, modeled via a 
QQ pair coupled to a constituent, non-relativistic string 
(a "flux-tube") (HIHIllll^l. It is therefore desirable to 
compare the "universal" decay picture of QQ and QQg 
mesons emerging from the Coulomb gauge, as illustrated 
in Fig. ^ with the exotic meson decay phenomenology 
based on models of flux tube breaking. 

In this paper we will investigate hybrid meson decays 
in a relativistic. Coulomb-gauge motivated description. 
This is a necessary first step towards understanding the 
phenomenology of normal meson decays. In the follow- 
ing section we discuss the construction of the Coulomb 
wave function for mesons and exotic hybrid mesons. For 
exotic hybrids we concentrate on the J^^ = 1 ^ quan- 
tum numbers which are expected for the lightest exotic 
multiplets [2^ |2^ [s^, |^ • In Section III we discuss the 
role of relativistic effects and give numerical predictions 
for various decay modes. In Section IV we consider de- 
cays of normal mesons, in particular the p and &i. In 
Section V we summarize our results and present future 
plans. 




FIG. 1: Quark flow diagram responsible for strong decays of 
QQg hybrids mesons (upper) and QQ normal mesons (lower). 
In the upper diagram the dashed lines represent the confining 
non-abelian Coulomb potential. The gluons connecting the 
Coulomb lines represent formation of the flux tube, e.g.the 
gluon string in the ground state. The overall state depicted is 
enclosed by the solid oval. Emission of a quark-antiquark pair 
from a gluon leads to a re-adjustment of the gluon strings and 
formation of the normal mesons in the final state. In the lower 
diagram the hybrid meson state appears as an intermediate 
state in a normal meson decay, which is assumed to proceed 
via mixing of a QQ pair with virtual excitation of a gluonic 
string and its subsequent decay. 



II. MESON AND HYBRID WAVE FUNCTIONS 

Recently lattice data has become available for static 
QQ potentials with excited gluonic flux 0, . In Born- 
Oppeinhemer approximation the lightest exotic hybrids, 
which by definition do not mix with the ground state QQ 
configuration, correspond to states built on top of the 
first excited adiabatic potential. The gluonic configura- 
tions can be classified according to symmetries of the QQ 
system, similar to the case of a diatomic molecule. The 
strong interaction is invariant under rotations around the 
QQ axis, reflection in a plane containing the two sources, 
and with respect to the product of parity and charge 
conjugation; thus each configuration can be labeled by 
the corresponding eigenvalues denoted by A — 0, 1, • • •, 
Y = ±1, and PC = ±1 respectively. In the ground 
state the gluonic flux tube has A = and in the flrst 
excited state it has one unit of spin, A = 1. Furthermore 
in the first excited configuration, lattice simulations find 
PC = -l for the gluon cloud 9_, J£,^] 

In the Coulomb gauge the quantum numbers of gluonic 
states can be associated with those of the extra transverse 
gluon in the presence of a static QQ source. This is 
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FIG. 2: Two-body (left) and three-body (right) potential 
between the transverse gluon and the static QQ sources. 
The dashed line represents the expectation value of the non- 
abelian Coulomb potential. 



because the transverse gluons are dressed [19|, |33 , and on 
average behave as constituents with effective mass ~ 
600 MeV,34. .35j . Thus low lying excited gluonic states 
are expected to have a small number of transverse gluons. 
The flux tube itself is expected to emerge from the strong 
coupling of transverse gluons to the Coulomb potential. 
The transverse gluon wave function can be obtained by 
diagonalizing the net quark-antiquark-gluon interactions 
shown in Fig.|21(in addition to the gluon kinetic energy). 

The three-body interaction shown in Fig. |21 is of par- 
ticular relevance. If the gluon is in a relative S'-wave with 
respect to the QQ system it has a net PC = +1, which 
is opposite to the lattice results for the first excited con- 
figuration. However, the transverse gluon has a gradient 
coupling to the Coulomb potential, thus a P-wave trans- 
verse gluon receives no energy shift from this coupling 
and the energy of the S'-wave gluon state is increased. 
In the Coulomb gauge picture, the shift of the S'-wave 
state via this 3-body interaction may be the cause of the 
inversion oi S — P levels seen on the lattice |3^ . A quan- 
titative analysis of this effect will be the subject of a 
separate investigation. In the following, when consider- 
ing the wave function of an exotic hybrid we will assume 
that the transverse gluon is in a P-wave relative to the 
QQ state. 



A. Mesons as QQ bound states 

In this analysis we do not solve the Coulomb gauge 
QCD Hamiltonian to obtain meson wave functions. In- 
stead, we use the general transformation properties un- 
der the remaining kincmatical symmetries (rotations and 
translations) to construct the states. Finally we em- 
ploy Lorentz transformations as for non-interacting con- 
stituents to obtain meson wave functions for finite center 
of mass (c.o.m.) momenta that are required in decay cal- 
culations. This is clearly an approximation which cannot 
be avoided without solving the dynamical equations for 
the boost generators jl^ . In the following section we 



examine some single particle observables to investigate 
this approximation. 

For a system of non-interacting particles the spin wave 
function is constructed as an element of an irreducible. 



unitary representation of the Poincare group [371, yM> 1^ 
I40I l4lj |. We will assume isospin symmetry m„^d = m 
and treat the light flavors first. The generalization to 
strange and heavier mesons will be given at the end of 
this section. 

In the rest frame of a quark-antiquark pair the quark 
momenta are given by 

= (P(m, q), q), 1^^ = (i?(m, -q), -q), (1) 

and the normalized spin-0 and spin-1 wave function cor- 
responding to J^'-'' = ^ and J^'~^ = 1 are given by 
the Clebsch-Gordan coefficients. 
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,(7g|0,0) 



(2) 



and 



V2 



(3) 



respectively. These can also be expressed in terms of 
Dirac spinors as 



1 



V2r, 



-u{q, CTg)7^w(-q, CTg) 



and 



(q: = 0,Crg,(Tg) 



2q' 



nio 



2m 



(4) 
(5) 



where mgg is the invariant mass of the QQ pair, rrigg = 
E{m, q) -I- P(m, — q) = 2\/m2 + q^, and e'(Agg) are the 
polarization vectors corresponding to spin 1 quantized 
along the z-axis. The boosted spin functions for mesons 
are summarized in Appendix A. The wave function of a 
QQ system moving with a total momentum Igg = Ig+lg 7^ 
is given by 



*g|' (q, Ag, Ag-) = J2 k - 0, CTq, CTq) 



xn*(i/2)/ 1 ^n(l/2)^_ 1 A 



(6) 



where the Wigner rotation matrix D^j^^^\q,P) corre- 
sponds to a boost with (3^ = P/M. One can show (see 
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Appendix A) that the corresponding wave functions can 
be written in terms of covariant amplitudes, 



1?97 ^q^ ^q) 



1 



V2r, 



-u{l„Xqh'vikAq) (7) 



and 



*<j|''(qi Iggj \^ ^q) 



£fj.O-qq, Xqq) 
V^rUgq 



u(^qAq) 



K - I'i 



Triqq + 2m 



vilq,Xq), 



(8) 



where e^^{lqq,Xqq) are obtained from (0, e*(Aqg)) through 
a Lorentz boost with Pj = Iqq/mqq. Obviously Ig — 
A(0 ^ lg,-)q and Iq = A(0 ^ lgg)(-q). 

By coupling Eqs. Q or ||SJ| for Igg = with one unit of 
the orbital angular momentum L — 1 and then making 
the boost in Eq. ©, one obtains respectively the spin 
wave function for the quark-antiquark pair with quantum 
numbers J^'^ = 1"' , 



qq 



V2mqq{lq,\^ 



-uilq,Xq)^'v{lq,Xq)Y^X,,{ci), (9) 



or for = 0++, 1++ and 2++, 



99 



-u{lq,Xq) 



-qq 



- 

''g ''9 



xw(lg, Ag)e^(lgg, A)Yi;(q)(l, A; 1, l\J, A,,), 



(10) 



with q = A(lgg 0)lg. In order to construct meson 
spin wave functions for higher orbital angular momenta 
L between the QQ pair one need only to replace Yn with 
Yli- For consistency, we should add the constant factor 
loo to wave functions with L ^ 0. 

Now we can proceed with the construction of me- 
son states characterized by momentum P, spin Xqq, and 
isospin /. The tt (/ = 1) and 77 {/ = 0) states ( J^"^ ^ 
^) are constructed in terms of the annihilation and 
creation operators: 



|M(P,/,/3)) - E / 



d^Pqd^ Pq 



all A,c,/ 

'N{P) V3 V2 



^2E{m,Pq)2E{m,Pq) 

*99(P9' P9' \' >^q)i^L{mqq{pq,Pq)/fl) 



2{E{mq,Pq) + E{mq, Pq)) ■ [2^)'' 5' [Pq + Pq - P) 



qqWq, I'qi ^Vg, ysq j ^y_Lyn .qgy^jq , ^jq j , fA. , "p^A,/,c,"p5A,-/,-C5 l^)' 



di 



(11) 



In the above 5'gg is the spin-0 wave function of Eq. 
written explicitly in terms of the momenta Pg and Pg 
instead of the relativistic relative momentum q and the 
c.o.m. momentum P — Igg, given by 



(q P)P 



E{m, q) 



E{mqq, P){mqq + E{mqq, P)) 



and 



(12) 



(q P)P 



E{m,q) 



E{mqq,P){mqq + E{mqq, P)) 



(13) 

In Eq. l(TT)) . c, / and I3 denote respectively the color, fla- 
vor, and third component of isospin. ip^ represents the 
orbital wave function resulting from the quark-antiquark 
interaction that leads to a bound state (meson). We as- 
sume that this function depends only on the invariant 
mass mgg. The normalization constant N (with P = |P|) 
is fixed by 

(P,a|P',a') = {2n)^2E{mM,P)S^{P --P')Sac.', (14) 

where ruM is the meson mass and a represents spin and 
isospin. The parameter is a scalar function of the me- 



son quantum numbers. Finally, F{I, I^) is a 2 x 2 isospin 
matrix (/ = 1 for u and f — 2 for d): 



F(0,0) = /, F{l,h)^a'e\h). 



(15) 



The flavor structure of the ry state (as well as other 
isospin zero mesons) was chosen as a linear combination 
-^{\uu) + \dd)), although in general those states are lin- 
ear combinations cos(0)[|uu) -I- |dd)]/'\/2-l-sin((/))|ss). The 
\ss) does not contribute to the amplitude of the decay of 
the TTi and therefore may be neglected in calculations, 
provided this amplitude is multiplied by a factor cos(0). 

Similarly the p{I = 1) and lo, 4){I = 0) states ( J^*^ = 
1 ) are given by Eq. (|ll(l . but instead of '^qq in Eq. 

O one must use *g|' given in Eq. (jSJ). The 61 (/ = 1) 

and hi (/ = 0) states (J^*^ — 1^ ) contain the wave 
function of Eq. Q. Finally, the a (/ = 1) and f {I = 0) 
states (J^^ = 0, 1, 2++) correspond to Eq. 

The above results can be straightforwardly generalized 
to the case where ruq and rriq are different, for example 
to decays into mesons with one strange quark (/ = 1/2). 
The spin wave function for a quark-antiquark pair in a 
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= state is 

(Ig, Ig, Aq, \q) = U{mq, Iq, \q)^^v{mq, Iq, Xq) 



■-qq 



(16) 

where fhqq = ^Jra^q — {ruq — rriq)^, and the if-meson 

states are given by Eq. with an appropriate defi- 

nition of the matrix F. In Eq. lfTT|l . if "^qq in Eq. lfTB|l is 

replaced by 



V2r. 



-u{mq,lq,Xq 



i-qq 



- l'^ 
'g 'g 



Xw(TOq,lg, Aq)e^(lqq, Aqg), 



(17) 



then one obtains the iir*-meson states (J^ = 1^). The 
wave functions of strange mesons with non-zero angu- 
lar momentum (such as = 1"*") can be obtained by 
coupling with the spherical harmonics. 

So far we have treated mesons as non-interacting QQ 
pairs. The interaction between a quark and an antiquark 
enters through the Hamiltonian H = and the boost 
generators of the Poincare group M"*. It is possible to 
produce models of interactions for a fixed number of con- 
stituents that preserve the Poincare algebra following the 
prescription of Bakamjian and Thomas |22.5 32J. Unfortu- 
nately such a construction does not guarantee that phys- 
ical observables, 6.(7. current matrix elements and decay 
amplitudes will obey relativistic covariance. In any case 
one deals with phenomenological models of the quark dy- 
namics, therefore we follow the common practice of em- 
ploying a simple (Gaussian) parametrization of the or- 
bital wave functions with one scale parameter /i related 
to the size of the meson, 



V'L(m,,-^) = e-<^/«'^'. 



gluon with momentum +Q, the total spin wave function 
of the hybrid is obtained by coupling the QQ spin-1 wave 
function of Eq. © and the gluon wave function ( J^*^ = 
1 — ) to total spin S' = 0, 1, 2 and J^^ = 0++, 1++, 2++ 
states respectively. The exotic meson wave function with 
J^'-' = 1 ^ is then obtained by adding one unit of orbital 
angular momentum between the gluon and the QQ pair: 

KT9iS)iKh,^a)^ E *gg-'(q,-Q,Ag,A,) 

xYu{Q){l, Xqq; 1, a\S, M)Di'>*{Q){S, M; 1, l\l, Xe.). 

(19) 

The spin-1 rotation matrix D^^\ representing the gluon 
spin wave function, relates the transverse gluon states in 
the helicity basis a{= ±1) to the basis described by spin 
Xg quantized along a fixed z-axis. The Clebsch-Gordan 
coefficients and the spherical harmonic in Eq. 119|) can 
be expressed in terms of the polarization vectors, and the 
action of the rotation matrix on the gluon states results 
in replacing e*(Ag) with e*(Q,Ag), where 



eliQ, Xq) = {XgW^ ^Q^Q^). 



(20) 



(18) 



Using the construction of spin wave functions summa- 
rized in Appendix B, the corresponding normalized wave 
functions are then given by 

*gg-;(S) = E -Q' ^'?)C(S)(Q, \q, Ag, A,,), 

(21) 

where the spin states C( s) El- 1211 are given by 



C(s=o) - V^[e*(Agg) • Xq)][Q ■ e(Ae.)], (22) 



In the non-relativistic limit (where the Wigner rotation 
may be ignored) all the spin wave functions for regular 
mesons simply reduce to the Clebsch-Gordan coefficients 
that we started from. 



C(S=1) 



-[e*(A,,-)xe:(Q,A<,)]-[Qxe(Ae.)], (23) 



B. Hybrid mesons as QQg bound states 



27 
104^ 



Q.[e*(Agg)®e:(Q,A3)]-e(Ae.), (24) 



As we discussed previously, in the exotic hybrid meson 
wave function the gluon is expected to have one unit 
of orbital angular momentum with respect to the QQ 
pair. In the rest frame of the 3-body system, where the 
QQ pair moves with momentum — Q and the transverse 



with Q = Q/IQI and {A(g)B)i-i = 2A(,Bj) - |(5ij(A-B) 



It is easy to show C(2) — 3/\/13(C(i) — fC(o))- The loss of 
linear independence is directly related to the transversity 
ofe^ 

The hybrid state in its rest frame is given by 



kl(/3, Ae^)) 



E 

all \,c,f 



1 



(pT£)q d^Pqd^Q 

{2Tr)^2Eq2EqEg 



i27r)^2iEq + Eq- + Eg)SHpq 



Q) 



'f,f. 



V2 



X 



*gg-s (Pg, Pg-, Ag, A,-, Ag)V'i bl, fdl 



,|0), 



(25) 



6 



k 


-p+1 




q-1 




-k 


-iP-q 



matrices. 

The orbital angular momentum wave function for the 
QQg system should depend only on the invariant masses 
TOgg, TTiqqg, and we will again introduce an exponential 
function 

(26) 

In the non-relativistic limit only that part of the tti spin 
wave function corresponding to the QQ pair is reduced 
to a Clebsch-Gordan coefficient, whereas the functions C, 
given in Eqs. H22I - I 24|l remain unchanged. 



FIG. 3: Coulomb gauge description of normal meson decays. 
The QQ component of the meson mixes with the QQg state 
with a subsequent decay of the transverse gluon into a QQ 
pair. The vertical solid line represents the sum over all QQg 
intermediate states. 



where the spin wave function '^qqg was given in Eq. H21|) 
for 5 = 0, 1, 2, and the orbital wave function V'^ depends 
only on niqq and the invariant mass of the three-body 
system, niqqg 



Here denotes the effective mass of the 



gluon and Eg = E{mg, Q), while Ac^c, are the Gell-Mann 



Non-exotic hybrids or QQg components of 
normal mesons 



As discussed in Section 1, in Coulomb gauge the decay 
of a normal QQ meson is expected to proceed via its 
QQg component with the gluon dissociating into a QQ 
pair, as shown in Fig. |21 The QQg component of the 
wave function is obtained by integrating the QQ state 
over the amplitude of transverse gluon emission from the 
Coulomb line 0|, shown by the vertical dashed line in 
Fig. 121 which gives, 



(2^)3 



^2ujg{q-\)AE 



r 



In Eq. (|?7|l V{p, q) is given by a product of the two 
Faddeev-Popov operators (represented by the dashed line 
in Fig. ^ modified by a Coulomb kernel vertex correc- 
tion. The transverse gluon couples to the Coulomb line 
via a derivative coupling which produces the momentum 
dependence of the numerator. The denominator is given 
by the difference between the energy of the QQ state and 
the energy of the QQg hybrid (non-exotic) state in the 
absence of mixing between the two j42| | . 

The spin-orbit structure of the QQg component can 
be inferred from Eq. H27(l . The momentum vector in the 
numerator can be coupled with the spin-orbit compo- 
nent of the \E'gg wave function and later coupled with the 
J^*-^ = 1 transverse gluon. For example for the p me- 
son the QQg component can be expanded in a basis of 
the ao, «!, 02 - like QQ wave functions all having spin one 
and one unit of orbital angular momentum between the 
quark and antiquark, coupled with the transverse gluon 
wave function ec(q — 1, Ag) to give J^'^ = 1 . Specifi- 
cally for the p meson one obtains, 



E 



(Q: — "Qj Ag, Ag) 



(28) 



Here ^'g^^'" are the Gq, ai and a2 QQ wave functions for 
J = 0,1,2 respectively, and the gluon helicity a = ±1. 
The normalized wave functions for p, similar to those of 
TTi in Eq. (|21|l . are then given by 



= E \> Ag)C(J)(Q, A, Ag, Ap), 



^ qqa{J) 
where 

C(J=0) 

C(j=i) 



— [e*(A)-q][e:(Q,Ag).e(A 



(29) 



(30) 



J=2) 



327r 
27 
leOTT 



[e*(A)xq].K(Q,Ag)xe(Ap)],(31) 



e:(Q,Ag).[e*(A)®q]-e(Ap). (32) 



Here q denotes again the quark momentum in the rest 
frame of the QQ pair. The most general wave function 
will be given by a linear combination of the three compo- 
nents listed above, and it can be calculated from Eq. l(77)>. 
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We will also study decays of the bi meson which are of- 
ten used as a testing ground for decay models. The QQg 
wave function with J^^ — 1^ , 1=1 quantum numbers 
requires the QQ to have tt or tt2 quantum numbers. The 
corresponding total wave functions are given by 



x^i?;(Q)(>/,A,,;l,(7|l,A,,) 



with ^qq"" being the tt {112) QQ wave function for J = 
(J = 2). Thus the normalized spin wave functions are 
given by 



(33) 



J2 -Q' Ag)C(J)(Q, A, Ag, Ap), 



where 



and 



C(J=0) 



[6:(Q,Ag).e(AbJ] 



(35) 



In the constituent basis used here the single-particle 
quark and antiquark orbitals correspond to states of 
massive particles with relativistic dispersion relations, in 
which a running quark mass is approximated by a con- 
stant constituent mass, 



=22] (2.)32i^(m,k) ["('^'^)'-^-^ + 



(38) 



Similarly the gluon field A'^' is expanded in a basis of 
transverse quasi-gluons, with a single particle wave func- 
tion characterizing a state of mass irig ^ 600 MeV, 



(34) y: I 



[e,(k,A)a^«;,+e:(-k,A)aL^^Je' 



(27r)32i;(TOg,k) 



(39) 

In Eq. H37() . g is the strong coupling constant, later chosen 
to be of the order 10, corresponding to as = 0{\). The 
decay matrix element 



C(J=2 



27 



^q-K(Q,A3)®6(AfcJ].q. (36) 



III. EXOTIC MESON DECAYS 



A. Relativistic 5*1 model 



(Mi(P),M2(-P)|i/|7ri) = (27r)353(p_p)A(P) (40) 

(where M denotes a meson) determines the decay ampli- 
tude A. 

For decays of tti into vrjy or tt^i, the spin part of the 
amplitude A is proportional to 



We assume that the transverse gluon in the tti creates 
a quark-antiquark pair and the hybrid decays into two 
mesons with momenta P and —P. Since the quark pair 
is emitted in the S = \, L = Q state this decay mechanism 
is also referred to as the ^S'l model. The Hamiltonian H 
of this process in Coulomb gauge is given by 



E 



d^^^,^f{-K){g-f ■ A='(x))V^e2/(x) 



A 



E ^M(5)(Q, A„ A,-, Xg,Xe.)ei{Q, Ag), (41) 



where 



(37) and 



Tr 



— m){p — m) {^^ -t- 



It' 



mgg(p, 1) -I- 2m 
I 



(43) 



In this expression p and I are respectively the four- 
momenta of the quark and the antiquark in the tti, 
whereas r and k are the four-momenta of the quark and 
the antiquark created from the gluon. If jir^ = /x^ then 
Aitrj — 0. The 1 ^ state is also found to have vanish- 
ing decay amplitude into two pions because of a relative 



minus sign from isospin that makes both terms cancel. 
Thus we find that the 1 ^ isovector does not decay to 
identical pseudoscalars. This is a relativistic generaliza- 
tion of a sy mmetry found in other non-relativistic decay 
models . Since tt and 77 are to a good approximation 
members of the same flavor multiplet, in the quark model 
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one typically finds their orbital wave functions to be simi- 
lar, i.e.,iJ,Tr ~ l^n, resulting in a small tti — > rjn decay rate. 
Of the two decay channels nr] and irbi , the latter will be 
favored. However, the parameters /x need not be similar 
for two mesons with different radial quantum numbers, 
making corresponding channels significant. For decays 
of TTi into np, irfj or rjaj (J = 0,2), the spin part is 



J 



proportional to 

E te(Q' \' ^fl' Ae.)e^,(Q, A,)e''*(-P, A), 

As 

(44) 

where A is the spin of the second meson, and 



{t> + m) (7" ^ ) (/ - m) (7" \ ) (/ + m)7^' - m)7' 

V mgg(p, 1) + 2m/ V mgg(r,l) + 2m/ 



(45) 



If /tip = /Xtt the amplitude of the decay into 7rp does not 
vanish (unlike tttj) and this channel can be favored. The 

same holds for tti tt/j and tti ^ ryaj. For the decay 
TTi — ^ puj, the spin part is proportional to 

D^""" E (Q' ^fl)^'c(Q, A,)e''* (P, Ap)e^* (-P, A„), 

(46) 



where 



mgg(p, k) + 2m 



)(i' + m)(n 



mqq{p, 1) + 2m 
I 



)(/-m)(i 



mgg(r, 1) + 2m 



^(/-|-m)7 

m 



If pp = fjLu then A = and the hybrid will not decay 
into p and w. Because both parameters p are expected 
to be of the same order, the channel tti — > pu will not be 
favored. 

Finally, for decays into strange mesons, one should 
use the above spin factors (depending on the quan- 
tum numbers), with a small modification resulting from 

= fc^ = m^. The amplitudes Akki {Sqq = 0, 1) behave 
similarly to A^^i and A^f^, whereas Akk' is like A^p. 
Therefore the former will be dominant and the latter is 



J 



expected to be much smaller. 

The width rate for a decay into a final state with orbital 
angular momentum L is equal to 



Po 



327r2 



mf 



-aliPo), 



(48) 



where Pq is defined by £'(mi, Pq) + E{m2, — Po) 
The partial wave amplitudes are given by 



Top 



aL{P) = E / ^(P,Ai,A2,Ae,)(Ji,Ai; J2,A2|J,A)yiM(P)(J,A;L,M|l,Ae,)df^, 



J,\,M 



(49) 



with dO. being the element of the solid angle in the di- 
rection of P, while m and A are respectively the masses 
and spins of the outgoing mesons. 



B. Non-relativistic limit 

The non-relativistic limit is obtained when the quark 
masses arc large compared to the quantities /i and Pq- 
This is equivalent to ignoring the Wigner rotation and 
taking non-relativistic phase space. In the orbital wave 



9 



functions, however, we must keep next to leading order 
terms that depend on momenta, otherwise the amphtude 
would become divergent. For tti — > ttt], nbi the dominant 
term has the form 



-32m* S'^ 



(50) 



while the other components are much smaller. Therefore 
the spin factor of Eq. (|41|l vanishes for S — I, and for 
5 = it tends to 8^/2l/TTm*Q'e\Xe^). From the above 



it follows F 



(S=l) 



and F 



TTi 7rp, 7r/i we have 



{S=2)/^{S=0) 



4/13. For 



(51) 



and the other components are much smaller. Therefore 
the spin factor of Eq. H44|) vanishes for S* = 0, whereas 



for S 

Thus F(5=o) 

for TTl 



1 it tends to — 8-\/6/7rim^(^ 
and F(s=2)/r(s=i) 
pLD we have 



eJ'(Ae,)e'=*(A)e''=^'. 
9/13. Finally, 



32to'*((5*^<5'=' - 5'^5^^ + 5'^5^^) (52) 



(the other components are again much smaller). There- 
fore the spin factor of Eq. 146|l tends for 5 = to 
8y/2A/^rri*Q^e\\ex)e^*{\p)ti*{\u,), and for S* = 1 to 
8V67^m4(Q^5^"'= - QJ5''=)e"(Ap)e^*(A<.)e'=(Aex). 

We can straightforwardly understand the difference in 
amplitudes coming from the spin wave function. If we 



assume ra,- 



''pi 



Hr, = Hp and TTibi = TO/j = m/2, ^fc^ 



M/i ~ M /2 (the second condition for masses is satisfied to 
a good approximation), then one obtains 



A — — A V — Ir 



and 



(53) 



(54) 



where At^^i, ^Trbi are taken for S* = and A-^p, A^^f-^ ^ for 
5' = 1. The relation between and A-^b-^^ (or between 
At^p and ^jr/i 2 ) is rnore complicated and depends on the 
orbital angular momentum wave functions ipL and ip'j^. 
If /ip — fi-^ then in the non-relativistic limit tti will not 
decay into np. Therefore the width for this process is 
expected to be much smaller than that of irbi, assuming 
the parameters /ip and /x^ are very close to one another. 
Analogous calculations can be made for the decays of 
TTl into strange mesons. If fix* — I^-k then in the non- 
relativistic limit Akk' — 



Numerical results 



Our model contains the following free parameters: the 
quark masses m; the effective gluon mass nig] the size 



FIG. 4: The pion weak decay constant as a function of 
the pion momentum p for m = 306 MeV and = 221 MeV. 



parameters /i for the wave functions; and the strong cou- 
pling constant g. The wave function parameters are con- 
strained by the pion decay constant and the elastic 
form factor f^. They are respectively given by 



and 



(0|^^''(0)K'=(p))=/rf^<5.fe, (55) 



(^''(p')l^^'^(0)|7r'=(p)) = F,{p^+p'^)ie,,k ■ (56) 



The axial current A'^''(0) is equal to ?/'c/(0)7^75(T*-!/'c//2 
and the vector current y^-' (0) to i^cfi'^h^o^^cf/^.. By 
virtue of the Lorentz invariance /tt is a constant, whereas 
is a function of Q-^ = — (p — p')^- As was discussed 
previously, it is not possible to construct the currents and 
wave functions with a fixed number of constituents in a 
Lorentz covariant manner. Thus the current matrix ele- 
ments are expected to violate Lorentz covariance. This 
will be reflected, for example in different values of 
obtained from spatial and time components of the axial 
current (rotational symmetry is not broken) . Even if we 
replaced the factor E{mM,P) in Eq. (|14|) by 1, it would 
be very difficult to find generators of the Poincare group 
that satisfy the exact commutation relations. Thus our 
model with exponential orbital wave functions will not 
be completely invariant. The resulting form factors will 
depend on the frame of reference and in order to obtain 
Ftt{Q^ = 0) = 1 one typically employs the time compo- 
nent /i = and works in the Breit frame. For other light 
unflavored mesons we do not expect much variation on 
the ground state wave function and we choose the scale 
parameter p to be of the same order as /i^. For kaons 
the same procedure fits and pK- 

We will assume that the mass difference between tt and 
p arises only from spin. Therefore we can write 



TOM = TOM + k{si ■ S2), 



(57) 



where M denotes either meson and mm is its 'averaged' 
mass. Substituting = 140 MeV and rrip = 770 MeV 
we obtain rfiM = 612 MeV, and thus for the constituent 
quark masses we choose = md = rhM/2 = 306 MeV. 
A similar relation can be used for the K and K* mesons, 
leading to ffiK = 792 MeV and mg = ffiK ~ to„ — 
486 MeV. The 'averaged' masses should be used in the 
normalization constants. 
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Mil 



FIG. 5: The pion electromagnetic form factor as a func- 
tion of the momentum transfer for m = 306 MeV and 
^r: = 221 MeV. The experimental values are from Ref. [i^ . 



TreZ nrel ) 


S = Q 


S = 1 


5 = 2 


7r6i(1235) 


S 


150(259) 


<1(0) 


44(80) 


D 


<1« 1) 


<1(0) 


< 1« 1) 


7r/i(1285) 


S 


<1(0) 


20(33) 


14(23) 


D 


<1(0) 


<1«1) 


< 1« 1) 


71-/2(1270) 


D 


<1(0) 


<1(<1) 


< 1« 1) 


7rp(770) 


P 


3(0) 


<1(0) 


1(0) 


7fA:*(892) 


P 


1(0) 


<1(0) 


<1(0) 



TABLE I: Relativistic (non-relativistic) widths in MeV for 
decays of the 7ri(1600). 



The weak decay constants can be used to fit the pa- 
rameters and ^K- Because of the Lorentz covari- 
ance breaking mentioned before, they become a func- 
tion of the meson momentum and we choose them to 
be equal to their experimental values at rest. Thus set- 
ting /^(O) = 93 MeV and /^(O) = 113 MeV, leads to 
^JL^, = 221 MeV and ^ik = 275 MeV. The momentum 
dependence of in our model is presented in Fig. 01 
which shows ^ 20% difference between the value of /^r 
calculated for the meson at rest and for the meson with 
momentum approaching the light cone. In Fig. |5l we 
present F^{Q^) calculated with the wave function param- 
eters obtained above and compared with data. There is 
good agreement for small momentum transfer; the dis- 
crepancy for larger indicates the lack of sufficient 
high momentum components in our wave function. The 
strong coupling constant at this scale is approximately 
= 10, and we take nig — 500 MeV for the gluon ef- 
fective mass. Around this value the Coulomb interaction 
between quark and antiquark appears to be linear. 

Next we proceed to discuss relativistic effects in tti ex- 
otic meson decays. In Tables I and II we present the rates 
for various decay channels. The numbers in parentheses 
correspond to calculations using the non-relativistic for- 
mulae, where S denotes the total spin of the QQg compo- 
nent of the TTi wave function. For all unflavored mesons 
and the tti , the value of the parameter fi was taken equal 
to /Itt, and for all strange mesons /i was set equal to ^k- 
This assumption makes the widths for the channels Trry, 
TTjy' and puj identically equal to zero. 

In Fig. we compare relativistic and non-relativistic 
predictions for the width for the decay tti irbi as 
a function of the mass of the light quark m. It also 
shows the semi-relativistic values, which include rela- 





5 = 


S = 1 


5" = 2 


7r6i(1235) 


s 


48(70) 


< 1(0) 


13(22) 


D 


1(2) 


< 1(0) 


< 1(< 1) 


7r/i(1285) 


S 


<1(0) 


7(11) 


5(8) 


D 


<1(0) 


2(< 1) 


1« 1) 


7r/2(1270) 


D 


<1(0) 


2(< 1) 


1« 1) 




P 


2(0) 


< 1(0) 


< 1(0) 


»7ai(1260) 


S 


<1(0) 


13(22) 


9(16) 


D 


<1(0) 


1« 1) 


< 1(< 1) 


rya2(1320) 


D 


<1(0) 


1« 1) 


< 1(< 1) 


i^A'i(1400) 


S 


127(45) 


< 1(0) 


39(14) 


D 


<1« 1) 


< 1(0) 


< 1(< 1) 


KKi{1270) 


S 


<1(0) 


11(4) 


7(3) 


D 


<1(0) 


< 1(< 1) 


< 1(< 1) 


KK* (892) 


P 


1(0) 


<1(0) 


<1(0) 



TABLE II: Relativistic (non-relativistic) widths in MeV for 
decays of the 7ri(2000). 




FIG. 6: Relativistic, non-relativistic and semi-relativistic 
widths for tti irfei in the S-wave state as a function of 
m, for S = and m^^ = 1600 MeV. 



tivistic phase space and orbital wave functions, but no 
Wigner rotations. The ratios of non-relativistic to rel- 
ativistic (and semi-relativistic to relativistic) values are 
shown in Fig. [7| 

From these results it is clear that fully relativistic re- 
sults are significantly different from non-relativistic ones. 
There are two sources of this difference: the Wigner rota- 
tion which introduces relativistic coupling between spin 
and spatial degrees of freedom in the wave functions, and 




FIG. 7: Ratios of non-relativistic to relativistic, and semi- 
relativistic to relativistic width rates for tti —> nbi in the 
S-wave as a function of m, for 5 = and m^x ~ 1600 MeV. 
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different relations between energy, momentum and the 
invariant masses (in the phase space and orbital wave 
functions). For realistic quark masses, both corrections 
appear to introduce corrections as large as 10% and thus 
should be included in phenomenological models. 



IV. NORMAL MESON DECAYS 

In this section we will calculate the widths of the de- 
cays p —> 27r and bi — > ttuj. These are the dominant 
decay channels (accounting for almost 100% of the total 
width) and their values are well known from experiment, 
and therefore can be used to test the model presented in 
this work. 

As discussed previously, in Coulomb gauge the decay 
of a normal meson is expected to proceed via QQ mixing 
with the QQg hybrid component followed by gluon dis- 
sociation to a QQ pair. However, the common approach 
to normal meson decays is based on the ^Pq model where 
QQ pair creation is described by an effective operator 
that creates the pair from the vacuum, in the presence of 
the normal QQ component of the decaying meson. We 
will first discuss the role of relativistic effects in the "^Pq 
model and then compare with predictions based on the 
Coulomb gauge picture. 



A. The decay p tttt 
We start from the '^Po Hamiltonian 



H = AJ2J d'x4/(x)Vc/(x), 



(58) 



where A is a mass scale that can be fixed by the absolute 
decay width, and is expected to be of the order of the 
average quark momentum. In the decay matrix element 
the spin factor is proportional to 



r 



x(/-m)(/-m) e\\p), 



(59) 



with the same notation used in the previous section. In 
the non-relativistic limit the above expression tends to 
32to^(p' — P*)e'(Ap). The expression for this process is 



much simpler than for decays of tti because there is no 
gluon. 

The decay of a p treated as a gluonic bound state 
has a similar structure to the decay tti — > ■nrj, and pro- 
ceeds via the interaction given by Eq. The 
spin factor is the same as in Eq. H41|l with S — J 
and Xex replaced by Ap, but the functions ^ are now 
given by Eqs. (30-32). In the non-relativistic limit this 



factor will become 
J=0, -8^9/27rm'^q'e^{Xp){5'^ + 



'{XpW^ - Q'Q^) for 
Qi) for J=l. and 
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\2-k\ P II 59(195) 1 7(15) 1 12(45) 1 58(75)] 

TABLE IIL Relativistic (non-relativistic) widths in MeV for 
the decay p(770) tttt for A = /x^. 



-8^27/107rm4gV^(Ap)i(7(5*^' - Q^Q^) for J=2. None of 
these functions vanishes, but only two of them remain 
linearly independent. 

In Table III, we present numerical predictions for the 
widths of this decay. The experimental value of the width 
for p — > 27r is 149 MeV. This number can be used to fit 
the free parameter A, the coupling constant g, or the 
hybrid scale parameter Hex'- 



B. The decay b\ —f nu 

This process is a better test for the different decay 
schemes, because the ratio of D-wave to S-wave widths 
is independent of the values of A and g. We will begin 
with the picture of the bi as a QQ bound state and the 
decay Hamiltonian of Eq. H58|) . i.e., the ^Pq model. The 
spin factor is proportional to 



Tr 



( / -|- m) {/k — m) — m) ( / — m) 



)]e;(-P,A.), 



(60) 



with the same notation as for p 2tt. In the non- 
relativistic limit this expression tends to 32m^{p^ — 
P')e"(A^). The spin factor for the decay of a bi treated 
as a gluonic bound state is given by Eq. I|41|) in which 



B''^ = Tr 



{ Jk — m){ p — m){ jl — m) 



mgg(r,l) -I- 2m 



(/ -I- to)7"' 



(61) 



r 



with S = J, Xex replaced by Ab^, and the functions ^ 
given by Eqs. (35-36). 



The numerical predictions for the decay widths are pre- 
sented in Table IV. The experimental value of the total 
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7ra;(782) 


S 


16(15) 


52(82) 


D 


17(42) 


<1« 1) 



TABLE IV: Relativistic (non-relativistic) widths in MeV of 
the decay fei(1235) 7rcj(782) for A = fi.^. 



width for the process bi ttcj is 142 MeV, and the ex- 
perimental ratio of the D-wave to S-wave widths is 0.08. 
Our predictions give a value less than 0.02 for this ratio 
in the model, and close to 1 for the ^Pq decay. The 
real mechanism for this decay appears to lie somewhere 
between the two predictions. The QQg wave function 
component of the bi wave functions used here is that of 
Eq. (|34|l . corresponding to a QQ pair with pion quan- 
tum numbers. For a QQ pair with 7r2 quantum numbers 
the numerical value for the width is much smaller than 
1 MeV for the S-wave and approximately 1 MeV for the 
D-wave. The ratio D/S is respectively 230. In the non- 
relativistic limit one obtains similar results. It is clear 
that treating bi as 7r2 +g increases dramatically the D/S 
ratio; therefore this may be an important component of 
the wave function. The relatively small values of the 
decay widths of a bi with 7r2 quantum numbers (L—2) 
compared to those of a bi with pion quantum numbers 
(L=0) remind the situation for the process tti — + nbi, 
whose D-wave width was small compared to the S-wave. 

V. SUMMARY AND OUTLOOK 

In this work we studied relativistic effects for the de- 
cays of normal and exotic mesons, and discussed a new 
picture of meson decays in the Coulomb gauge point of 
view. In Coulomb gauge the gluon degrees of freedom 
are physical. Since they carry color, isolated gluons do 
not appear in the physical spectrum, and colorless exci- 
tations of QQg states are expected to be suppressed by 
a mass gap of the order of 1 GeV. This is what lattice 
QCD studies find for the energy of gluonic excitations 
in the presence of QQ sources. Since strong decays are 
expected to proceed via gluon decay into a QQ pair, the 
Coulomb gauge provides a natural framework for disen- 
tangling the dynamics of bound state formation (via a 
static Coulomb potential) and the decay of the gluonic 
component of a state. 

This work led to two important conclusions. First, 
numerical results showed significant relativistic correc- 
tions arising from spin-orbit correlations introduced by 
Wigner rotation. The widths calculated using fully rel- 
ativistic formulae are in general larger than the corre- 
sponding values calculated with no Wigner rotation (by 
a factor of roughly 10%), but smaller than completely 
non-relativistic rates. Some decays that are suppressed 
in the non-relativistic limit (for example tti 7rp assum- 



ing identical orbital wave functions for tt and p) acquire 
non-zero amplitudes in the relativistic case. 

The second conclusion from this work is that the 
lightest exotic meson, the tti, prefers to decay into two 
mesons, one of which has no orbital angular momentum 
while the other has L = 1 (the so-called S + P selec- 
tion rule) [2^. Thus this selection r ule, also found in 
other models seems to be quite robust Q . Some decays 
{ttt], puj, KK*) are suppressed by symmetries in orbital 
wave functions, or the assumption that the parameter 
p should be almost equal for mesons with the same ra- 
dial quantum numbers. We have also noticed that, for 
decays where two waves are possible, the rates for the 
higher partial waves are larger in the ^Pq than in the 
^Si model. However, there is one caveat that needs to 
be explored further. There are several components to 
the QQg normal meson wave functions and if there are 
sizable contributions from wave functions with large rel- 
ative angular momentum in the QQ system, it is possible 
to obtain large amplitudes for high partial waves. This 
has in particular been illustrated in the case of the D- 
wave/S-wave ratio for the bi ttw decay. 

The process tti irp seems to be suppressed and this 
agrees with the S + P selection rule |4J|. However, some 
models predict larger values for its width . It is possi- 
ble that these are increased by the final state interactions 
between the outgoing mesons. The most important con- 
tribution may come from the process nbi — > vrp, that 
proceeds through an u exchange. Calculations of this 
contribution will be the subject of future work. 
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VI. APPENDIX 

A. Boosted spin wave functions for mesons 

The polarization vectors corresponding to spin 1 quan- 
tized along the z-axis are given by 

The Wigner rotation matrix corresponding to a boost 
with = P/Af is given by 
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{E{m, q) + m){E{M, P) + A-/) + P • q + icr • (P x q) 



^2{E{m, q) + to) {E{M, P) + M){E{m, q}E{M, P) + P • q + mM) 

I 



AA' 



whereas 5(1^^ 0) is the Dirac representation of the 
boost taking 1^ to q and 1^ to — q: 



Sik, ^ 0) 



I'qq, '■qq 

• Igg 



E{mqq, \qq) + TOgg 



(63) 



(64) 



We make use of the relations: 



(65) 



^ <X'^(q, lgg)«^(q, <rq) = i.t(i^^ Ag)^t(i^_ ^ 0), 



and the formulae: 



(66) 



e.{\qq) = A;(0^ lgg)e,(A„-), 

one obtains Eq. ijSJ from Eq. jSJ and Eq. ©. One can 

derive the spin wave functions for other mesons in similar 
fashion. 



B. Spin wave function of the tti 

The transverse gluon states in the helicity basis a are 
related to the states in the spin basis Ag (quantized along 
a fixed z-axis) by 

\Q,Xq) =^Z?«;(<^,0,-0)|Q,a), (67) 



where 9 and cf) are respectively the polar angle and az- 
imuth of the gluon momentum direction. For the gluon 
polarization vector we obtain 

eciQ,Xq)= E <l*('/',^,-0)4(Q,^), (68) 

cr=±l 



where the helicity polarization vectors are given by 

Using the unitarity of the matrix D'^^ one can show 

4(Q,Ag)6r(Q,a) = Z?«;(Q), (70) 
and with the help of the identity 

er(Q,^)4(Q,^) - 5''~Q'Q' 

we arrive at the result: 



eliQ,\q) = e\XqW^ -Q^Q^^ 



(71) 



where = QV|Q|. 

The Clebsch-Gordan coefficients and the spherical har- 
monic in Eq. p9|) can be expressed in terms of the po- 
larization vectors H62|) . for example: 



(1,A';0,0|1,A) = e*(A')-e(A), 
(1,A';1,A|0,0) = e*(A')-e*(A), 



(72) 
(73) 



(1,A';1,A"|1,A) 
YuiQ) 

Therefore we obtain: 



^/2 



[6*(A')xe*(A")]-6(A), (74) 



-e{l) ■ Q. 



(75) 



^(1, Xqq; 1, Agio, 0)yiz(Q)(0, 0; 1, /|1, Ae.) = [e*(Agg) • e*(Ag)][Q • e(Ae.)], 



(76) 
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^(l,A,,-;l,A<,|l,s)yi/(Q)(l,5;l,;|l,Ae.) = [e*(A„-) x e*(A5)] • [Q x e(Ae.)], (77) 

l,s 

^(l,A,,-;l,Ag|2,s)yiz(Q)(2,s;l,;|l,Ae.) = Q • [e*(A,,-) e*(Ag)] • e(Ae.). (78) 

Ls 



The action of the rotation matrix I?^^^ on the gluon states to the spin wave functions given in Eqs. (|22I24I) . 
results in replacing e'(Ag) with e*(Q, Ag) and that leads 
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